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1 Introduction 
1.1 Frozen boundaries 

Given two polynomials P{z,w) and Q{z,w), one can study solutions of the 
system 

Q{e-''z,e-yw) = 0. ^ ^ 

as a function of a point {x,y) G M?. A solution z{x,y) and w{x,y) of (Hj) 
solves a first order quasilinear PDE 

^ + ^ = 1, P{z,w) = 0, 
z w 

which is closely related to the complex Burgers equation and arises in the 
theory of random surfaces, see [12]. The singularities of z{x,y) and w{x,y) 
occur when (^Q) has a multiple root, that is, when the two curves in (^Q) are 
tangent. In the random surface context, this marks the boundary between 
order (e.g. crystalline facet) and disorder, called frozen boundary. 

For example, let P and Q define rational curves in MP^ of degree 2 and 
4, respectively, positioned with respect to the coordinate axes of MP^ as 
illustrated in Figure ^ Points where P and Q intersect the axes may be fixed 
so that the frozen boundary will be inscribed in a hexagon as illustrated in 
Figure El The probabilistic meaning of this curve is the following. 

The well-known arctic circle theorem of Cohn, Larsen, and Propp j3] 
states that the limit shape of a typical 3D partition contained in a cube 
has a frozen boundary which is a circle. More precisely, projected in the 



1 



Figure 1: Curves P and Q for the frozen boundary in Figure El 




Figure 2: An example of frozen boundary 

(1, 1, l)-direction, the cube becomes a hexagon and the (1, 1, l)-projection 
of the frozen boundary is an inscribed circle. Now suppose we additionally 
weight each 3D partition vr by a weight which a product over all boxes 

in TT of some periodic function oi i — j and j — k. If the period equals 2, the 
frozen boundary will look as in Figure |21 The coefficients of the curve P are 
obtained from the periodic weights, while the coefficients of Q are fixed by 
boundary conditions. See |lH] for a detailed discussion of this procedure in 
general. 

The geometry and, especially, the singularities of frozen boundaries are 
of considerable interest. Note that these are curves of some complexity: in 
our current example frozen boundary has degree 28 and genus 27. 
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1.2 Log- front 



This motivates the following definition. Given two curves P,Q C. (C^)^, 
consider 

R= {{a,b)\ Ta,b ■ Q is tangent to P} C {Cf (2) 

where 

it the dilation {z, w) i— {az, bw) and tangency at a singular point means that 
for some branches of the two curves in question their unique tangent lines 
coincide. For simplicity, we assume that tangency occurs only at isolated 
points, that is, no component of P is a dilate of a component of Q. 

We call the curve R the log-front of P and Q. The reason for such name 
wiU be explained below. We will denote the construction Q by 

R = P^Q. 

The obvious property 

Q ^ P = R{a'\b-^) 

and (jH) below justify such notation. 

Clearly, R is an algebraic curve. Various ways to compute its equation 
R{a, b) = will be discussed in Section 13.11 The frozen boundary from 
Section [TTT] is given by i?(e^',e^) = 0. 

1.3 Symmetry between Q and R 

The definition of R can be recast into several equivalent forms. Let V denote 
the hypersurface 

V = {P{az,bw) = 0} G X (C^)2 

and let iTat and tt^w be the projections from (C^)^ x (C^)^ to the respective 
factors. 

The tangency in definition (j2)) can be rephrased by saying that R is formed 
by critical values of the map 

TT,,: S = Vn7r-^'{Q)^{C-f (3) 

of complex surfaces. In other words, the surface 

^at (R) n < (Q) = i? X Q 
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is tangent to the hypersurface V along their intersection, showing a certain 
symmetry between the roles of Q and R for fixed P. In particular 

gcp-(p-g). (4) 

The multiplicity, with which Q occurs in the right-hand side of (0} equals 
the degree of the logarithmic Gaufi map of P, see below. 

Yet another way to say the same thing is that R is the envelope of the 
family of curves 

{P(az,M=0}Q(,,^)^„, (5) 

indexed by points of Q 

1.4 Classical constructions 

Among examples of the operation (j^J there are the following two classical 
constructions. 

First, let P be a general line, for example, 

P(z, w) = z + w + 1 . 

Then its dilates Tab ■ P form an open set of the dual projective plane and 
hence R is an open set of the dual curve Q^, namely, its intersection with 
(C^)^. In this case, (jH) becomes an equality. 
Second, the additive analog of © 

P-Q= {(a, 6) \ Q{z - a,w -b) =0 is tangent to P} C (6) 

is a limit case of (j21). For analytic plane curves, constructions P ^ Q and 
® are, in fact, equivalent by taking the logarithms. When P is a circle of 
radius t 

P{z,w) = z^ + -t'^ 

the real locus of P — Q contains the front at time t of a wave that was 
emitted at time zero from all real points of Q and is propagating with unit 
velocity. Such curve is called a wave-front. It is this example that motivates 
the general term log-front. 

Pliicker formulas relate singularities of a curve Q to the singularities of its 
dual Q^. A formula of Klein further constraints the singularities of the real 
locus of Q^. Analogous formulas for wave fronts were obtained by O. Viro in 
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PB] . The goal of this note is prove an analog of Pliicker and Klein formulas 
in the general case. 

In probabilistic applications, the curve P is a real algebraic curve of a very 
special kind, namely, it is a Harnack curves. See and the Appendix for 
a discussion of the properties of Harnack curves and ^] for connections 
with probability. Harnack curves have many remarkable features that general 
real plane curves lack. As it turns out, the assumption that P is Harnack is 
also essential for our derivation of Klein-type formulas for R. 

It may be noted here that both our formulas and their proofs involve 
nothing but elementary geometry of plane curves and would have been, no 
doubt, obtained by Pliicker and Klein had their seen a need for them. 

1.5 Acknowledgments 

We are grateful to C. Faber, R. Kenyon and O. Viro for many useful discus- 
sions. Some of our formulas generalize unpublished results obtained jointly 
with C. Faber. 

2 Preliminaries 

There are many good books on plane algebraic curves, see e.g. |SII2H]- With 
the probability audience in mind, we collected in this section an explanation 
of some basic notions that will be used later. 

2.1 Newton polygons 

Let P C (C^)^ be a curve with the equation 



This is a refinement of the degree of P. The polygon Ap is defined only up to 
translation by a lattice vector once we treat P C (C^)^ as a geometric curve 
rather than a polynomial. Let |Ap| denote the number of lattice points in 




By definition, its Newton polygon is 



Ap = Convex Hull {(i, j) G Z'^\pij ^ O} . 



Ap. 
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The Newton polygon defines a projective toric surface 



r=\(z^wj).. ... I c (7) 

of degree 

degT = 2Area(Ap) . 

From now on, P will denote the closed curve in T defined by P{z, w) = 0. 
This is a hyperplane section avoiding the torus fixed points, which are the 
only possible singularities of T. We call the components T>e of 

edges E of Ap 

the boundary divisors and the points of 

dP = PndT 

the boundary points of P. For simplicity we assume throughout the paper 
that the boundary points of both P and Q are smooth. 

By construction, P intersects T>e in l^'l points counting multiplicity, 
where \E\ is the length of E in lattice units. The multiplicities of these 
intersection points define a partition X{E) of \E\ for every edge E of Ap. 
We will call a lattice polygon with such additional partition data a marked 
polygon. 

2.2 Amoebas 

The image Log(P) of a curve P under the map 

Log{{z,w)) = (log|2;|,log|u7|) 

is called the amoeba of P. The boundary points of P give rise to the so- 
called tentacles of the amoeba Log(P), see for example Figure El which shows 
logarithmic images Log(RP) of the real loci of two plane curves. 

The directions of the tentacles are the outward normals to the sides of 
the Newton polygon. For example, Newton polygons of curves from Figure 
El are plotted in Figure El The number of tentacles in a given direction is the 
number of points of P on the corresponding boundary divisor. 
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Figure 3: Newtons polygons of curves in Figure El 
2.3 Logarithmic Gaufi map 

Since (C^)^ is an abelian group, the tangent spaces at all points of (C^)^ are 
canonically identified with tangent space ^(11) at the identity (1, 1) G (C^)^. 
Mapping the tangents to P to their images in P(T(i 1)) gives the logarithmic 
Gaufi map 

where P is the normalization of P. In coordinates, 

7P((^,^)) = ^4^. 

w-^P 

aw 

Note that the values of jp at the boundary points of P are determined by 
the slopes of the corresponding edges of Ap. 

The degree of 7p can be computed as follows, cf. jTTll2n]. Recall that the 
multiplicity m{p) of a point p G P is the multiplicity with which P intersects 
a generic line through p. Another important characteristic of a singular point 
is its Milnor number fi{p), see j22]- It niay be defined as the local intersection 
number of -^P = and ^P = at p. 

az aw 

Proposition 2.1. We have 

deg7p = 2|A^| + \dP\ - 2 - + m(p) - 1) . (8) 

where |Ap| is the number of points in the interior of the Newton polygon Ap 
and \dP\ is the number of boundary points of P not counting multiplicity. 
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Clearly, only singular points of P contribute to the sum over p G P in (jSI) 
since for a nonsingular point p we have m{p) = 1 and fi{p) = 0. 

Proof. If P is smooth and transverse to the boundary then the proposition 
follows from Kouchnirenko's formula ^H] (also using Pick's formula for the 
area of a lattice polygon) . 

A singular point of P subtracts fi{p) + m{p) — 1, which is the same com- 
putation as for the ordinary GauB map, see e.g. Theorem 7.2.2. in [2H1- The 
effect of tangency to the boundary is determined similarly, see e.g. Section 
lOl below. □ 

The critical points of 7p are known as the logarithmic inflection points. 
Note that, for example, a cusp is not a logarithmic inflection point, but a 
sharp, or ramphoid, cusp (locally looking like = y^) is. By the Riemann- 
Hurwitz formula applied to jp, the count of logarithmic inflection points 
with multiplicity is 

^ logarithmic inflection points = 2deg7p — x{P)- (9) 
2.4 Geometric genus and adjunction 

The geometric genus of P is defined as the genus g{P) of its normalization 
P. It may be computed as follows, cf. |THj . 

9{P) = \^p\-\Y.iKp)+m-^) (10) 

where P{p) is the number of branches of P through p. This is known as the 
adjunction formula. 

The number |(/i(p) + /3(p) — 1) is a nonnegative integer which vanishes 
if p is a smooth point of P. It equals 1 for cusps and nodes, so if the curve 
has only those singularities, we get 

# cusps + # nodes= |A^| -^(P). (11) 

Combining © and ^ we get 

^(m(p) - f3{p)) = - deg7p - x{P \ dP) . (12) 

pGP 
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Remark 2.2. A geometrically-minded reader might appreciate the following 
alternative proof (or rather a rephrasing of the proof) of IT^ 

Suppose that P \ dP is smooth. We may deduce the equation 

deg^P = -x{P-^dP) (13) 

by the following application of the maximum principle, cf. P^. Let A : — > 
M be a linear map which we may choose with a generic slope of the kernel. 
The function A o Log is a pluriharmonic function on (C^)^ and thus restricts 
to a harmonic function h on deg7p = —x{P \ dP). Thus all critical points 
of h are of index 1. 

Furthermore, h exhibits P as a cobordism between P_ C dP and P+ C 
dP, where the splitting dP = P+UP_ corresponds to h^^{+oo) and h~^{—oo). 
Thus, the number of critical points of h equals x{P ~ P+i P-) = x{P ^ dP) 
(the latter equality follows from additivity of Euler characteristic). On the 
other hand the number equals deg7p, so we get (fT^. Each singularity p of 
P_ C dP subtracts ^(/i(P)-/?(p)- 1) from -x(P) and ^(^(P) -m(p) - 1) 
from deg7p. 

2.5 Nodal and cuspidal numbers of a plane curve 

The right-hand side of (jl2p may be interpreted as the total cuspidal number 
c{P) of P. It equals the number of cusps if P has only cusps and ordinary 
multiple points. An example of formula ()12j) may be seen in Figured In the 
bottom half of Figure 01 we see two inflection point come together to form a 
point of the form y = x^. For the dual curves, which have the same genus 
and the same degree of the logarithmic GauB map, we see a merger of cusps 
into a Eq singularity = which has the cuspidal number equal to 2. 



Figure 4: Two cusps coalescing in a Eq singularity. 
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Another useful number is the nodal number h{P). First we define the local 
nodal number bs{P) at every singular point s G P as the sum ^ ^{Bj,Bk) 

of the intersection numbers over all distinct pairs of branches Bj, Bk via s. 
Clearly, bs{P) = if and only if s is a locally irreducible singular point of P. 
If s is an ordinary node of P then bs{P) = 1. The nodal number b{P) is the 
sum of the local nodal numbers over all singular points of P. 

In many applications the numbers c{P) and b{P) play the role of the total 
numbers of cusps and nodes even if P has higher singularities. 

2.6 Integration w.r.t. Euler characteristic 

Because of the equation 

x{A nB) = x{A) + x{B) - x{A n B) 

the Euler characteristic may be viewed as a signed finite-additive measure. 
If / is a function on X taking finite many values, its integral with respect to 
the Euler characteristic is defined by 

/ fdx = J2y^(f''(y^y 

y 

Calculus of such integrals was developed by O. Viro |23- Applications of 
this calculus to the classical Pliicker and Klein formulas may also be found 
in gll. 

More generally, given a map 

F -.X 

one may push-forward F^f oi f by the formula 

[F.f]{y)= [ fdx. 

Under additional hypotheses on F and /, this operation has natural func- 
torial properties like (FG)^ = F^G^, see e.g. Section 7.3 in |28]. The case 
of main importance for us will be when F is a nonconstant map of smooth 
curves, in which case, 

F.l = degF-brF, 

where hi F is the branch divisor of F, that is, the sum of all critical points 
of F with their multiplicities. 
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3 The equation of R 



3.1 as a Resultant 



Consider the curve TZ defined by the equations 




0. 



(14) 



in (C^)^ X (C^)^. We have TZ G S, where the surface S was defined in 
In general, TZ may have several components. The Wronskian in (fT^ vanishes 
at any singular point of g of Q, hence such a point contributes a translate 
of P to TZ. The multiplicity of this component equals m(g) + — 1, see 
the proof of Proposition 12 .11 Symmetrically, singular points of P contribute 
copies of Q to TZ. The remaining components correspond to actual tangency. 
The curve R is the projection of these components to the (a, 6)-plane. 

The equation of R may be found by eliminating z and w from (jl4|) . fol- 
lowed by factoring off spurious components caused by singularities. Grobner 
basis algorithms give one way to perform this elimination. We find that in 
practice it is easier and faster to use resultants for this computation. 

Consider the polynomial 



The equation Ri{a, b, z) = cuts out the image Si of S under the projection 
along the z direction. This projection may create singularities, specifically 
Si may intersect itself along curves. Consider 



Its zero locus consists of critical points of the projection of Si along the z 
direction, together with contributions of singularities. Those can be identified 
and removed by factoring i?2 as they all occur with multiplicity greater than 
one. For example, the image of double point curves of Si will occur in R2 
with multiplicity two. 

Exact factorization is possible and rather effective for polynomials with 
rational coefficients. In probabilistic applications, the coefficients may be 
known only approximately, making exact factorization impossible. This is 
why it is useful to know the Newton polygon of R, which will be determined 
below. Also note that sometimes, albeit rarely, components of R may occur 
with multiplicity in R2 as, for example, in (0)). 



Ri{a,b, z) = resultant^t, (-P(az, few), If;)) . 
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3.2 Multiplicity of tangency 



Given two curves C and C and a point p define their tangency multiplicity 

at p as 

(C A CX = (C ■ CX - mc{p) mcip) 

where (C • C')p is the local intersection multiplicity of C and C at p and 
mc{p) is the multiplicity of the point p on C. Note that (C A C")p = if 
p ^ C or if C and C intersect transversely at p. Also note that the tangency 
multiplicity is additive over branches of both C and C . 
We have the following 

Proposition 3.1. The multiplicity of a point (a, b) & R is the total tangency 
multiplicity of P{az, bw) = and Q{z, w) = 0, assuming they are not tangent 
at infinity. 

In other words, the local cuspidal number at each branch of a singular 
point of R is equal to the corresponding tangency multiplicity minus one. 

Proof. Since multiplicity is a local quantity, we can work with the additive 
version (jH)). We may work with each branch of the singularity individually. 
Suppose P and Q are tangent ai z = w = and that w = is their common 
tangent. In this case the origin will belong to P — Q and 6 = will be the 
corresponding tangent. We need to compute the total multiplicity mji of all 
branches of R corresponding to this tangency. It equals the number of roots 
bi of R{e, b) = with bi ^ where e is a fixed small nonzero number. For 
b = bi the curves 

P{z + e,w + b) = 0, Q{z,w) = 0, (15) 

are tangent. 

Consider the intersection U of the curve defined by ()15p in the {z,w,b)- 
space with a neighborhood of the origin. We will compute the Euler charac- 
teristic of U in two different ways. Viewing f/ as a degree mp mq branched 
covering of the z line, we get 

xiU) = -mp mq + mp + mq , 

where mp and mq is the multiplicity of the origin on P and Q, respectively. 



12 



On the other hand, we may view [/ as a branched covering of the 6-hne. 
It has degree {P-Q)o, where {P-Q)o is the intersection multiphcity of P and 
Q at the origin. We have 

X{U) = {P ■ Q)o - mp{mQ - 1) - mQ{mp - 1) - niR 

because each tangency in (|15|) corresponds to a simple branchpoint. This 
concludes proof. □ 



4 Newton polygon of log-front 
4.1 Formula for Ar 

Let Ap and Aq denote the marked Newton polytopes of P and Q. We 
denote by — Aq the reflection of Aq about the origin, with the corresponding 
marking. 

We say that two edges E and F are opposite and write E ^ F, ii their 
outward normals point in the opposite direction. As we will see, in the 
simplest case when no edge of — Aq is opposite to an edge of Ap, the Newton 
polygon Aji is simply the Minkowski sum 

deg(7Q) ■ Ap + deg(7p) ■ (-Aq) . (16) 

Every pair of opposite edges i^^ C Ap and F C — Aq makes the polygon Ar 
smaller than ()16|). 

Given two partitions A and fi, we define 

i 

= ^min(Ai,/ij) , (18) 

id 

where A' denotes the conjugate partition. The equivalence of (fTTj) and flTH|l 
is an elementary combinatorial fact. 
Given an edge E, let 

E 




is the unit (non-oriented) lattice segment in the direction of E. Note that 
(similarly to Newton polygons) such lattice segments are defined only up to 
translations by Z^. 
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Theorem 1. The Newton polygon of the log-front R is the unique polygon 
Afi satisfying 

An+Y, (A(^), A(F)) = deg(7Q) • Ap + deg(7p) ■ (-Aq) , (19) 

where the summation is over all pairs E G Ap and F C —Aq of opposite 
edges and X{E) are the corresponding partition markings, i.e. multiplicities 
of intersection with the boundary. 

Minkowski "subtraction" of segments, implicit in formula (fT^ . simply 
reduces the length of all edges in that direction by the given amount. 

4.2 Proof 

By definition, a boundary point of -R is a tangency of the curves P and Q, 
dilated with respect to one another by an infinite amount. In logarithmic 
coordinates, infinite dilation becomes an infinite shift. Shifted by an infinite 
amount, the amoeba Log(P) of P is either a union of parallel lines (tentacles) 
or empty, in which case there will not be any boundary points of R in this 
direction. This is illustrated in Figure 



Figure 5: Tangency to a tentacle creates a boundary point of a log-front 

Tangency to a tentacle occurs at points mapped by the logarithmic GauB 
map to the slope of the tentacle. Therefore, if Log(P) and Log{Q) do not 
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have tentacles in the same direction (that is, Ap and Aq have no parallel 
edges) then each tentacle of Log(P) contributes deg7Q points to dR. The 
multiplicity of the resulting point in dR equals the multiplicity of the corre- 
sponding point of dP. 

The effect of a tentacle-to-tentacle tangency may be studied in a local 
model, for example, 

In this case, 

k I 

7p(z) ~ -, 7q(z) -, z^l, 

z — 1 z — 1 

Equating 7p(az) = 7q(-2) gives a parametrization of R such that 

a~l, boc{z~lf-^, z^l. (20) 

This means that tentacles of P and Q pointing in the opposite direction 
(which means that k and I have opposite sign) produce a tentacle of com- 
bined multiplicity \k\ + |/| pointing in the same direction as the P-tentacle. 
Tentacles pointing in the same direction, by contrast, produce a tentacle 
of multiplicity \k — l\ pointing in one of the two directions. This rule can 
be phrased more naturally in terms of the reflected Newton polygon — Aq 
because the reflection flips the direction of tentacles. 

The polygons Ap and — Aq can, in total, have as many as 4 edges with 
the same slope. Let E, E C Ap and F, F C —Aq be such a 4-tuple of edges 
and let X{E), . . . denote the partition marking. Assume that the outward 
normals to E, F point in the same direction. By our computations, the 
polygon Ap will have an edge in the same direction as E and F of length 

|E|(deg7Q - iiKF)) - i{\(F))) + |F|(deg7P - i{X{E)) - i{X(E))) 

+ 5^(A(E). + A(F),)+ (21) 

id 

J]max(A(E)i - A(F)j, 0) + ^max(A(F)i - A(E),-,0) . 

id id 

Note that the middle line in ()2H1 cancels with a part of the first line. A 
further cancellation is obtained from 

-\E\ i{X(F)) + J2 max(A(E), - A(F)„ 0) = - J] min(A(E)„ X(F),) 
id id 
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using formula (fTH|). This concludes proof. 

Counting the boundary points of R without respect to the corresponding 
direction n^; gives us the following corollary. 

Corollary 4.1. The cardinality of \dR\ equals to 



deg7p|9AQ| + 7q|5Ap| - (A(E), A(F)). 



4.3 Examples 

Let P and Q be generic with Newton polygons from Figure El In this case, 
the Newton polygon A/j is the larger of the two polygons plotted in Figure IHl 
If P develops a tangency to the boundary then A^ shrinks and becomes the 
smaller polygon in Figure El What happens in this case, is that R becomes 
reducible with one component being a boundary divisor. The boundary 
divisor corresponds to tangency with P at the newly developed point of 
tangency to infinity. 




Figure 6: Possible A^ for Ap and Aq from Figure 



As another example, consider the case when P and Q are generic curves 
of degree dp and dq. In this case, both Ap and Aq are triangles with sides 
of length dp and dq, while Ap is a hexagon with sides of length 

dpdqidg - 1), dqdpidp - 1), ... , 
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cyclically repeated. In particular, when dp = 1 we get a triangle with side 
dqi^dq — 1), reproducing the very classical formula for the degree of dual 
curve. 

In the example in Figure El the degrees of logarithmic Gauss maps of P 
and Q are 4 and 10 respectively, hence Ar is a hexagon with sides 8 and 12. 

4.4 Reconstructing the curve Q by the log-front R 

It is instructive to have a closer look at the computation done in fl20j) in the 
case k = I. In this case, the point (a, b) does not escape at infinity, but it is 
still a pole of the logarithmic Gauss map of R (see also Section I^TD below) . 
In other words, the curve R is tangent to a = 1 at the corresponding point. 

In probabilistic applications, the curve P is given and one knows that the 
frozen boundary is compact and tangent to given lines. This allows to fix 
the real boundary points of Q and their multiplicities (which have to match 
the corresponding multiplicities for P). In particular, this is how the curve 
Q in Figure Q is determined from the requirement that the frozen boundary 
in Figure |21 is inscribed in a hexagon. 



5 Pliicker-type formulas for R 



5.1 Logarithmic Gaufi map of R 

Let R be the product of the normalizations P and Q over their logarithmic 
GauB maps: 

R (22) 




In English, a point of R corresponds to a tangency of a branch of P to a 
translate of a branch of Q. Note that -R^is a j)artial normalization of the 
plane curve R: we have the factorization R R ^ R for the normalization 
R ^ R. We claim that the logarithmic GauB map of R factors through the 
natural map R ^ R. 
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Proposition 5.1. The logarithmic Gaufi map for R is the composition 

R^R^¥^ 

Proof. The essential geometric content of tliis result is already implicit in 
(^. In coordinates, the claim is elementary to check on the open dense set 
where tangency is nondegenerate. It suffices to analyze the Gauss map of 
the additive analog P — Q of the log front P ^ Q. Let P and Q be given by 

y = P{x), y = Q{x). 

Then P — Q is parametrized by 

(x-s,p(x) -g(s)) 

where s = s{x) is a solution of P'{x) = Q'{s). The Gauss map of the above 
curve is clearly P'{x) = Q'{s). □ 

Corollary 1. 

deg -fR = deg 7p deg 7Q . (23) 
5.2 Geometric genus of R 

The geometric genus of R, or equivalently, the Euler characteristic x{R) can 
be computed by Riemann-Hurwitz formula applied to the map 7/j in ()22|) . 
To do this we need to be able to compare x(-R) and x(-R)- 

The map 7/? is ramified over the branchpoints of either 7p or 7q. Suppose 
that at a point p G P the map jp is given by 

x^x", z/ = l,2,..., (24) 

in a suitable local coordinate x E C centered at p. In other words, suppose 
that p is a logarithmic inflection point of multiplicity u — 1. Note that such 
point may be singular or non-singular point of P. This integer u will be 
denoted by i^(p). 

A point of i? is a pair {p,q) E P x Q such that 

7p(p) = iQiq) = xeF\ 

Let Xp and Xg be local coordinates at p and q as in (j24|) . The local equation 
of R is 

^'^(P) _ rf.'^ig) (rr ^ ^ (P2 
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This has gcd(i/(p), i^{q)) branches and hence produces gcd(z/(p), z/(g)) points 

'"i, • • • , rgcd{u{p),u(q)) e R 

that are mapped to x by the logarithmic GauB map. Each is a logarithmic 
inflection point of multiplicity lcm(z/(p), z/(g)) — 1. Note that while all G -R 
are logarithmic inflection points of the same multiplicity, the multiplicities 
of the corresponding branches of R may be different and are not determined 
by the numbers z/(p) and i^iq). 
From definitions 

X{R)= [ ldx= [ iipUiiQUdx. (25) 

We conclude 

Theorem 2. The Euler characteristic x(-R) is given by 

xiR) = xiR) + Yl (gcci(z/(p), z/(g)) - 1) 

7p(p)=7Q(g) 

where x{R) = ((7p)*1, (7q)*1)x = -2deg7p deg7Q+x(-P) deg7Q+x(Q) deg7p. 

The above inner product with respect to the Euler characteristic is defined 
as in ()25|). Note that generically the the branchpoints of 7p are disjoint from 
branchpoints of jq. In this case, the above formulas may be simplified as 
follows. 

Corollary 2. If tangency does not occur at two logarithmic inflection points 
then 

xiR) = X{R) = -2 deg 7p deg 7q + x{P) deg 7q + x{Q) deg 7p . (26) 
5.3 Nodes and cusps of R 

The total cuspidal number of R may now be determined using the formula 

Theorem 3. We have the following expressions for the cuspidal number cr 
of the resultant curve R = R{P, Q). 

c{R) = -x{R) + \dR\ - {x{P) + c{P) - \dP\){x{Q)+c{Q) - \dQ\). 
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Proof. This formula can be obtained as a straightforward combination of 
(|12|) and Corollary ^ Nevertheless it is instructive to prove it by the Viro 
calculus (see Section I2.6p to prepare a way for the real counterpart in the 
next section. 

Consider the family of translates Ts{P) parameterized by s G (C^)^. De- 
note with X{s) the Euler characteristic of the space of pairs {tp, tq) E P xQ, 
such that Tgitp) = tq. We have 

j X{s)dx{s) = {x{P)-\dP\){x{Q)-\dQ\). (27) 

{CX)2 

It may be viewed as a corollary of the Fubini theorem since each pair of 
points from P \ dP and Q \ dQ will appear in X{s) once. On the other 
hand we have 

j X{s)dx{s) = -x{R) + \dR\ 

(CX)2 

- c(P)(x(Q) - \dQ\) - c{Q){x{P) - \dP\) - c{P)c{Q) - c{R). (28) 

Indeed, we have x((C^)^) = while X{s) is constant for generic s E T 
(namely, by Bernstein-Kouchnirenko formula it is equal to Vo1(Ap,Aq)). 
However the value X{s) drops ii s E R and it drops further (by m — 1) if 
there is a singular point of multiplicity m at a branch of \ dR through s or 
if Ts{P) and Q have a tangency of higher order (which in turn corresponds to 
a cusp of R). We get the theorem as the combination of (P7|) and (|^. □ 

If P and Q are generic, the only singularities of R will be cusps and nodes 
and the nodal number of R may be recovered from the adjunction formula. 

For example, let P and Q be generic curves of degree dp and dq. In this 
case ^ 

X{R) = -dpdq{Adpdq - 3dp - 3dq) . 

Note that this grows very fast with dp and dq. Already when P and Q are 
generic conies, the geometric genus of i? is 9. 

The Newton polygon of R and its boundary were determined above in 
Section 14.31 Generically, R will not have multiple point on the boundary. 
For the number of cusps Theorem 13 produces 

# cusps = Sdj:,dQ — Gdpdq . 
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Accordingly for the number of nodes we get: 

# nodes = -y^ + 2rff,rf^ + -y^ - Sd^d^ - Sdj^d^ - 2rf^rf^ + 9dpdQ . 

For dp = 1, these speciahze to the classical Pliicker formulas. The last two 
formulas were obtained in (Zj. 



6 Klein-type formula for the log-front R 

6.1 Refinement of the nodal and cuspidal numbers 

For this section it will be important that both P and Q are defined over the 
field M of real numbers. In this case, clearly, the curve R is defined over M 
as well. 

An algebraic curve P C T is defined over M if and only if it is invariant 
with respect to the involution of complex conjugation conj : T T. The 
fixed point locus of this involution is the real toric surface RT C T. The 
real locus RP coincides with the intersection P fl RT. 

For a real curve P we may refine both the cuspidal number c(P) and 
the nodal number b{P) as follows. Let c^'^{P) be the number of cusps of P 
(counted with multiplicity as in Section that are real, i.e. contained in 
RTa. In other words, to get c^'^(P) we add over (real and imaginary) branches 
of real singular points of P the multiplicities of these branches diminished by 
1. We set c^'^(P) = c(P) - c^^(P). 

Let b^'^{P) be the sum over all singular points s of P in MT of the number 
of the pairs of conjugate imaginary branches of P. If s G RP is a singular 
point we may introduce b^^{P, s) as the sum of the local intersection numbers 
over all possible pairs of the real branches of P at s. Let b^'^{P) be the sum of 
b^'^{P, s) over all real singular points of P and let 6^™(P) = b{P) — 6^'^(P) — 
b^^{P). 

Thus we get the refinements 

c(P) = c^^{P) + c^^(P) 

and 

b{P) = bl'iP) + b^'{P) + b^'^iP). 

In the case when the only singularities of P are ordinary cusps and nodes 
the numbers c^^{P) and c^™(P) are the numbers of real and complex cusps 
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respectively while the numbers b^'^{P),b^'^{P) and 6^™(P) are the numbers 
of real hyperbolic (A^), real elliptic (Af) and imaginary nodes. 

It is convenient to define the boundary nodal number b{dP) to be equal to 
the sum of the multiplicities of the boundary points of P minus the simple 
cardinality of dP. (Note that since dP is 0-dimensional this boundary nodal 
number also works as the boundary counterpart of the cuspidal number.) 
Again we have the refinement b{dP) = b^'^{dP) + b^"^{dP) that counts real 
and imaginary boundary points separately. Thus b^^{dP) is the measure of 
nontransversality of the real locus MP to the boundary divisor of the toric 
surface T. 

6.2 Computations for R 

In general, the set of refined nodal and cuspidal numbers b^'^, b^'^, c^"^, c^™ for 
the log-front R is not determined by the corresponding sets for P and Q. 
Let us recall a well-known example illustrating this in the case of classical 
projective duality, i.e. if P is given by the polynomial z + w + 1: two real 
infiection points may disappear together with a bitangent real line. 

Example 6.1. Let MQi, MQa C MP^ be real quartic curves pictured in Figure 
[7| The first curve can be constructed by perturbation of the union of two 
ellipses while the second one can be constructed by perturbation of the union 
of four lines. We have b^^iRi) = 0, c^^(Pi) = 8, bf{R2) = 4, c^%Ri) = 0. 




Figure 7: Real quartic curves MQi and 'RQ2 whose projective duals have 
distinct b^'^ and c^'^. 

The classical Klein's formula JEj allows us to compute 2b^'^{R) + c^^{R) 
in the case when P is a line. To generalize this statement for a larger class 
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of real curves P let us look at the argument map Arg : (C^)^ x 
defined by 

ATg{z,w) = (arg(z), arg(w)). 
The image Arg(yl) G x is called the coamoeba or the alga of A C (C^)^, 

cf. ESI, m 

The real 2-torus x S*^ is a group which has Z2 x Z2 as its subgroup. 
Let T = S^x 5VZ2 X Z2 be the quotient group and P : x T = x 

X Z2 be the projection map. Note that T is itself a group isomorphic 
to X and the zero in this group is Ot = 1). 

We need to compactify the map (5 o Arg : P \ dP — > T. Note that this 
map does not extend to P. Let us consider P to be the real blow-up of P at 
the finite collection of points dP C P. Naturally, P is a closed non-orientable 
surface whose Euler characteristic coincides with that of P \ dP (in the case 
of non-compact spaces we use Euler characteristic for homology with closed 
support) Furthermore, we have a natural extension 

Algap : P^T 

such that Algap Ip^gp = (3 o Arg \p^Qp. 
Consider the following integral 




X{t)dxwT{t) 



t6Rr\Vcrt(IRr) 

(note that here we use the family Tf(P) only for t G WT). This integral is 
the real counterpart of the integral (j27|) from the proof of Theorem |21 where 
we defined the number X{t), but now we use only the real translations. 

As in Theorem E] the integral can be computed in two different ways. 

For this computation to depend only on visible characteristics of the 
curves P and Q we need to assume that P is a Harnack curve (see Appendix 
A). 

Let 

e(Ap,AQ)= Area(A'p + Ay, 

where Ap (resp. Ag) run over all possible sides of the polygon Ap (resp. 
Aq) and Ap + Aq stands for the parallelogram obtained as the Minkowski 
sum of the intervals Ap and Ag. 
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Proposition 6.2. If P is a Harnack curve and Q is any curve defined over 
R then 

/m = 2 Area(Ap)(x(Q) - 2bf{Q) - 6«^(9Q) - |9Aq|) + 

{2bf{P) + 6«^(9P))(26^^(Q) + b^'^idQ)) - e(Ap, Aq), 

Proof. Note that if two points sp E P and sq E Q are different by a trans- 
lation by t G WT then Algap(sp) = AlgaQ(sQ). Conversely, if Algap(sp) = 
AlgaQ(sQ) and either sp or sq is not in the exceptional divisor of the blowup 
map then there exists t G WT \ Vert(MT) such that Tt{P) and Q have sq as 
one of their intersection point. The Euler characteristic of the space formed 
by the pairs (sp,sq) with Algap(sp) = AlgaQ(sQ) and such that sp and 
Sq are from the exceptional divisors of the blowups P P and Q Q 
respectively is e(Ap, Aq). 

We have x((Algap)~^(a;)) = 2 Area(Ap) for any a G T \ {0^} by Lemma 
m Note that x((AlgaQ)-i(r \ {Ot})) = x{Q) - 2bf{Q) - b'^^dQ) = 
xiQ) - |<9Aq| - 2bf{Q) - b^%dQ). The contribution of a = Ot to the left- 
hand side of Proposition O is {2bf'{P) + b^%dP)){2bl%Q) + b^%dQ)) since 
MP is a closed 1-manifold and thus x(]RP) = 0. □ 

Let n be the number of vertices of the polygon Apg = Aq — Ap. 

Proposition 6.3. If P is a Harnack curve and Q is any curve defined over 
M immersed near the boundary o/MT then 

/m = Vol(Ap,AQ)(4-2n) 

- 2 Area(Ap)ci-(Q) - (26^'=(P) + 6^^(aP))c^^(Q) - {2bf{R) + CRe(P)). 

Proof. The Euler characteristic of M.T \ Vert(MT) is 4 — 2n. For a generic 
t we have X{t) = Vol(Ap, Aq) by the Bernstein-Kouchnirenko formula [2], 
|18j . This number gets decreased if Tt{P) and Q are tangent or if one of 
their intersection point is singular for a branch of Q (note that P is an 
immersed smooth curve since it is Harnack). The latter case contributes 
-2 Area(Ap)c^"^(Q) - (26f (P) + b^%dP))c^%Q). 

If Tt{P) and Q are tangent at a point in a non-real point then we have a 
bitangency since both P and Q are invariant with respect to the involution 
of complex conjugation. This contributes —2b^'^{R). The tangencies at real 
points contribute x(I^-R) = plus CRe(-R) where the tangencies of Tt{P) and 
Q are of higher order. □ 
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6.3 Klein's formula 



Combining Proposition 16.21 and 16.31 we get the following theorem. 

Theorem 4. If P is a simple Harnack curve and Q is any curve defined 
over M immersed near the boundary o/MT then 

2bf{R) + c^%R) = (4 - 2n) Vol(Ap, Aq) + 

2 Area(Ap)(|aAQ| - x{Q) + 2bf{Q) + b'^'idQ) - c'^^m 
- {2bf{P) + 6^^(9P))(26^^(g) + 6^^(9g) + c^^(Q)) + e(Ap, Aq). 

Corollary 6.4. Let MP, MQ e MP^ be two curves of degree dp and dg (re- 
spectively) not passing via (0 : : 1), (0 : 1 : 0), (1 : : 0). Suppose that P is 
a simple Harnack curve and Q is a smooth curve. Then 

2bf{R) + c^"(P) = 44 - 2dpdQ. 

Proof. We have n = 6, Area(Ap) = x{Q) = '^dq — dg, |i9Aq| = Srfg and 
e{Ap,AQ) = QdpdQ. □ 

Let us deduce the Klein formula ^7] in its classical form from Theorem 
m in the case when P is a line z + w + 1 = and is a curve of degree 
d in ]RP2 not passing via (0 : : 1), (0 : 1 : 0), (1 : : 0). In this case Ag 
is a triangle with vertices (0,0), (rf, 0) and {0,d). The line P is a simple 
Harnack curve of degree 1 with ap = 1 and thus we may apply Theorem EJ 
The polygon Ap is a hexagon in this case, so n = 6, while Vol(Ap, Aq) = d 
and e(Ap, Aq) = 6d. 

Note that the classical Klein formula computes 2b^'^{R) + c^'^{R) for the 
dual curve in MP^ while our formula does it for R C MT. The toric surface 
MT is the result of the blowup of MP^ at three points and, in general, such 
blowup might change the characteristics b^'^ and c^*^. Let us assume that Q 
intersects the boundary divisor CP^ \ (C^)^ transversely, so that this blowup 
is disjoint from Q. Then xiQ) = xiQ) ~ S"^- 

We get 

2bf{R) + c^{R) = -8d + 2bf{Q) - xiQ) + 3d - c^"'{Q) + 6d 

= d-c'^{Q) + 2bf{Q)-x{Q). 

Note that c^™(Q) = cq — c^^{Q) while the classical Pliicker formula [23] for Q 
implies that x{Q) + cq = 2d-d*, where d* is the degree of P C (C)^ C CP^ 
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(called the class of Q in Klein's paper [IHl)- Thus we get 2h^{R) +(^'^{R) = 
-d + d* + 2bf{Q) + c^%Q) and, therefore, 



which is the Klein formula in its original form 
6.4 Example 

Let us go back to the log- front from Figure |21 By Theorem 0] the sum of the 
real cusps of the log-front and twice the number of real solitary nodes equals 
24._ Indeed, we have n = 6, Vo1(Ap,Aq) = 8, Area(Ap) = 4, |(9Aq| = 12, 
xiQ) = 2, 2bf{Q) = h+{dQ) = c(Q) = while e(Ap, Aq) = 48. 

There are 6 real cusps visible on Figure |21 Furthermore, Figure |H1 shows 
images of the remaining 3 quadrants under the map (x, y) i-^ (log log \y\). 
There are 12 more real cusps. Thus by Theorem E] our log-front has 3 real 
solitary points. 



d - 2bf{Q) - c^"(Q) =d*- 2b] 



(/?)-c^^(i?). 



(29) 




Figure 8: Images in the other three quadrants 
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Appendix A: Harnack curves and their algae. 



In his 1876 paper [9^ A. Harnack produced for each d examples of algebraic 
curves of degree d in MP^ with ("^"^K*^"^) _|_ \ yqsI (topological) components. 
Furthermore, in the same paper he has shown that -)- i is the upper 

bound for the number of components of any curve of degree d in MP^ 

Among the curves constructed by Harnack there were some "canonical" 
curves whose topological arrangement in MP^ is especially easy to describe. 
Note that any component of a smooth curve is either contractible (i.e. bounds 
a disk in MP^) or is isotopic to MP^ C MP^. A contractible component is called 
an oval while the disk bounded by it is called the interior of the oval. The 
oval whose interior is disjoint from other ovals is called empty. 

For an odd d = 2/c + 1 there is a smooth algebraic curve of degree d that 

consists of k{2k — 1) empty ovals and a non-contractible component. For 

an even d = 2k there is a smooth algebraic curve of degree d that consists 

of {k — l){2k — 1) empty ovals and one other oval whose interior contains 
(fc-iKfc-2) _ ^^^^j^ _ g^p^y 

Let P C T be a curve given by a real polynomial with the Newton 
Polygon Ap and T be the toric surface corresponding to Ap. Since P is 
defined over M there is the real locus MP = P fl (M^)^ which coincides with 
the fixed point set of the involution of complex conjugation on P. 

Definition 6.5. A curve P is called a Harnack curve if for every (x, y) G 
the set Log~^(a;, y) O P consists of no more than two points. 

Remark 6.6. Earlier we called such curves simple Harnack curves to distin- 
guish them from other curves in the Harnack construction. However, by now 
we have convinced ourselves that these simple curves are the most beautiful 
in the Harnack series of constructions. We propose to drop "simple" from 
their name and call them Harnack curves. 

Harnack curves (from Definition I6.5|l exist for any convex lattice poly- 
gon A, cf. [TU]. In [inj it was shown that the topological type of the triad 
(MT; RP, dM.T) depends only on A if MP is a smooth Harnack curve trans- 
verse to infinity. 

In the case when A is a triangle with vertices (0,0), (rf, 0) and (0, rf) we 
have MT = MP^ and MP C MP^ is a curve consisting of ('^-^K'^-^) empty 
oval and one other component which is non-contractible if d is odd and an 
oval containing empty ovals if d is even. In |12, it was shown that 
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all such curve form a contractible subspace in the space of all real curves of 
degree d. 

Recently it was discovered that Harnack curves possess many extremal 
characteristic properties, among them are the following. 

• (EH) We have 

Area(Log(P)) = 7r^Area(A), 

if M.P is a Harnack curve with the Newton polygon A. In the same 
time by [21] for any curve A C (C^)^ given by a (not necessarily real) 
polynomial with the Newton polygon A we have 

Area(Log(A)) < 7r^Area(A). 

Furthermore, if Area(Log(y4)) = 7r^Area(A) then A can be translated 
(by a multiplication with some (2:0,^0) G (C^)^) to a Harnack curve, 
see EH- 

• (Cni) The curve Log(MP) C is embedded, does not have inflection 
points and contains #(Int(A)nZ^) compact components (called ovals) 
if MP is a Harnack curve. Each such oval comes from one of the four 
quadrants in (M^)^. The number of the ovals coming from the four 
quadrant equals to the number of lattice points (j, k) E I? in the 
interior Int(A) with given residue mod 2: j = jo (niod 2), k = fco 
(mod 2) (clearly there are four possible pairs of residues). 

Furthermore any curve with that many ovals of Log(MP) ^ P \ dP is 
Harnack if the remaining component of MP intersects the infinity 9MT 
in a maximal way (see P^). 

The only singularities of a Harnack curve MP are isolated double points 
(the singularities of type Af according to P) in (M^)^, see [21 . Note that 
even though MP has to be smooth near the boundary divisor dM.T it does 
not have to be transverse to the boundary divisor. An example of a singular 
simple Harnack curve of degree 6 in MP^ which is not transverse to the 
boundary divisor is sketched in Figure El 

The goal of this appendix is to give yet another characteristic property 
of the Harnack curves in terms of their algae. Suppose that P is a Harnack 
curve Recall that in the previous section we denoted with P the result of the 
real blowing up of P at the points dP and with 

Algap : P^T 
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Figure 9: A singular simple Harnack curve in MP^ and its position with 
respect to the boundary divisor. 

such that Algap Ip^gp = (3 o Arg |p^gp. Recall that Pq = Algap^(OT). 

Lemma Al. If P is a Harnack curve then the restriction of the map Algap 
to P \ Po 

P \ Po ^ T \ {Ot} 
is an unhranched covering of degree 2ylrea(Ap). 

Proof. The critical points of the map Arg|p are the points such that the 
image of the Logarithmic Gauss map 7p is in MP^ (see [T^). By [T^ 
7~^(]RP-^) = MP. The map Algap Ip^p^ is proper since Algap^(OT) = Pq. 

Since x(T \ {Ot}) = —1 we have the degree of the covering Algap Ip^p^ 
equal to — \ Pq)- Let us start with a generic P and study how does 
x{P \ Pq) changes when we deform P in the class of Harnack curve. Since P 
is a simple Harnack curve we have PxPq non-singular. The only singularities 
of P are the real isolated double points. Each such point contributes 2 to 
the Euler characteristic of the normalization of P, but also 2 gets subtracted 
when we remove this point from P. 

Note that P may also have the "boundary" singularity. This means that P 
is not transversal to the boundary divisor. In this case each point of tangency 
of order m between P and OT gives m — 1 points in P \ Pq. Thus, while 
x(P) gets increased by m — 1 in the case of such tangency (in comparison 
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with x{P) ill the transversal case) in turn x{P \ Po) g^ts decreased by m — 1. 
Thus in both cases the net effect of possible singularities on x(P\Po) is zero, 
so for the computation of the degree of our covering we may assume that P 
is smooth and transversal to the boundary divisor of T. In this case 

X{P X Po) = X{P) = X{P) = Area A. 
The last equality is a corollary of Khovanskii's formula □ 

We may compactify the set-up of Lemma Q to get the following Theorem 
describing the alga of a Harnack curve. Denote the blow-up of T centered in 
Ot with B :f ^T. 

Note that T is a surface with a natural involution induced by conj: 
(a,/5) 1-^ (— a, — /9), here we think of G 5*^ as arguments of complex 
numbers. Clearly, the 1-dimensional part of the fixed-point set of this invo- 
lution is the exceptional divisor of the blow-up. Denote it with M.T. 

For the compactifying theorem we need to blow up P even further. Recall 
that P is a smooth (real) surface equipped with an involution coming from 
complex conjugation. Note that Ot is the regular value of a map P\RP T 
and thus Pq C P is a disjoint union of embedded circles corresponding to the 
ovals of MP and some isolated points. Let n be the result of (real) blow-up 
of P at Pq. Clearly, a blow-up at a smooth submanifold of codimension 1 
does not change the surface P thus only blowups at the isolated points of Pq 
matter. Such points come either from isolated double points of P or from 
tangency of P with dT. Note that 

x(n) = x(P\Po). 

Theorem Al. There exists a map p completing the commutative diagram 



p 




The map p is a covering of degree 2 Area(Ap) which is equivariant with re- 
spect to the involution of complex conjugation defined on U. 

Furthermore, any real curve P such that Algap lifts to a covering n — >• T 
is a Harnack curve. 
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Proof. The point of WT is specified by a tangent line to T at 0^- Tlie 
logarithmic Gauss map 7p takes real values at M.P. Thus the tangent line 
at s G MP is real and gives (after multiplication by i) a tangent direction at 
Ot G T. We define the value p{s) to be this direction. At points of Pq \ RP 
the map is defined by the blowup itself (recall that Or is the regular value of 
the map P \ MP ^T). 

For the converse we note that by ^H] the amoeba map Log | p does not 
have critical points outside of MP since Algap does not have any. Thus 7p 
cannot have any critical points on MP. 

□ 
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